function b_value = compute_bvalue()
b1 = @(x,y) 1 - 2.*x;
b2 = @(x,y) 2.*x + 2.*y -1;
b3 = @(x,y) 1 - 2.*y;
gauss_points = [0.1012865073235, 0.1012865073253;
                0.7974269853531, 0.1012865073253;
                0.1012865073235, 0.7974269853531;
                0.4701420641051, 0.0597158717898;
                0.4701420641051, 0.4701420641051;
                0.0597158717898, 0.4701420641051;
                0.3333333333333, 0.3333333333333];
b_value = zeros(3,7);
for i = 1 : 7
    b_value(1,i) = b1(gauss_points(i,1),gauss_points(i,2));
    b_value(2,i) = b2(gauss_points(i,1),gauss_points(i,2));
    b_value(3,i) = b3(gauss_points(i,1),gauss_points(i,2));
end
end